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Abstract 

The three-dimensional Schrodinger's equation is analyzed with the help of the 
correspondence principle between classical and quantum-mechanical quantities. Sep- 



| aration is performed after reduction of the original equation to the form of the 

classical Hamilton-Jacobi equation. Each one-dimensional equation obtained after 
separation is solved by the conventional WKB method. Quasiclassical solution of 
the angular equation results in the integral of motion M 2 = (I + \) 2 h 2 and the exis- 
' fence of nontrivial solution for the angular quantum number I = 0. Generalization 

of the WKB method for multi-turning-point problems is given. Exact eigenvalues 
, for solvable and some "insoluble" spherically symmetric potentials are obtained. 

<^> ' Quasiclassical eigenfunctions are written in terms of elementary functions in the 

form of a standing wave. 



Q\ ' 1. Introduction 

Basic equation of quantum mechanics, the Schrodinger's wave equation, is usually 
solved in terms of special functions or numerically; for several potentials, the equation 
is solved exactly The general approach to solve the Schrodinger's equation for the 
qh solvable potentials^ is to reduce this equation to the equation for hypergeometric function 

or some special function. To do that one needs to find first a special transformation for the 
wave function and its arguments to reduce the original equation to the hypergeometric 
form. After that using certain requirements (defined by boundary conditions) to the 
hypergeometric function one can write the corresponding solution for the problem under 
consideration, i.e. the eigenfunctions and the corresponding eigenvalues. 

This is rather a mathematical approach to solve the eigenvalue problem in quantum 
mechanics; the corresponding methods and solutions of the wave equation for some po- 
tentials have been developed long before the creation of quantum mechanics. There are 
several features of this exact mathematical method that should be clarified from the phys- 
ical point of view. One of them is related to the S'-wave state. The radial Schrodinger's 
equation has no the centrifugal term for the orbital quantum number 1 = 0. From the 
physical point of view, this means that the problem does not have the left turning point. 
In order for the physical system to have a stable bound state (discrete spectrum) two 
turning points are required (see Ref. ffl, for instance). However, solving the radial wave 



"solvable" potentials we mean those models for which the eigenvalue problem can be reduced to 
a hypergeometric function by a suitable transformation. 
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equation for I > 0, one obtains energy eigenvalues for all / including / = 0. Another fea- 
ture is related to the angular dependence. The angular eigenfunction for the ground state, 
Y Q0 (6,ip) = const, i.e. no nontrivial solution exists. Meanwhile, as in the case of radial 
dependence, it might be a function with no zeroes of the type of a standing half-wave. 

There is another approach to the eigenvalue problem in quantum mechanics. This is 
the quasiclassical method which is well known and widely used mainly as the Wentzel- 
Kramers-Brillouin (WKB) approximation @-@] applicable in the case when the de Broglie 
wavelength, A = h/p (h = 2irh), is changing slowly. In several cases of interest, the WKB 
method yields the exact energy levels, however, its correct application results in the exact 
energy eigenvalues for all known solvable potentials. 

The quasiclassical method is based on the correspondence principle between classical 
functions and operators of quantum mechanics. The correspondence principle is used to 
derive the wave equation in quantum mechanics. In Ref. |J this principle has been used 
to derive the semiclassical wave equation appropriate in the quasiclassical region. It was 
shown that the standard WKB method (to leading order in h) is the appropriate method 
to solve this equation. 

In this work we solve the multi-dimensional Schrodinger's equation by the quasiclassi- 
cal method. Unlike known approaches, instead of modification of one- dimensional equa- 
tions obtained after separation, we analyze an original multi-dimensional Schrodinger's 
equation and reduce it to the equation in canonical form (without first derivatives). Sep- 
aration is performed after reduction of the equation obtained to the form of the classical 
Hamilton- Jacobi equation. We show that the main question of the exactness of the qua- 
siclassical method is tightly connected with the correspondence principle, i.e. the form 
of the generalized moments obtained after separation of the wave equation; the moments 
obtained after separation have to coincide with the corresponding classic generalized mo- 
ments. The quantization condition is written with help of the argument's principle in the 
complex plane that allows us to generalize the quasiclassical method for multi-turning- 
point problems and obtain the exact energy eigenvalues for all known solvable potentials 
and, also, for some "insoluble" problems with more then two turning points. 

The quasiclassical method reproduces not only the exact energy spectrum for known 
potentials but has new important features. One of the consequences of the quasiclassical 
solution of the mult i- dimensional Schrodinger's equation is the existence of a nontrivial 
angular solution at I — 0, Y^ KB (9, <p), which describes the quantum fluctuations of the 
angular moment. This method allows us to show apparently the contribution of quantum 
fluctuations of the angular momentum into the energy of the ground state. 

2. Exactness of the WKB method 

It is well known that the exact eigenvalues can be defined with the help of the asymp- 
totic solution, i.e. the exact solution and its asymptote correspond to the same exact 
eigenvalue of the problem under consideration. The asymptotic solutions in quantum 
mechanics can be obtained by the WKB method. Therefore the quasiclassical method 
should reproduce the exact energy spectrum. 

Intriguing results have been obtained with the help of the super symmetric WKB 
method (SWKB) 0-0, which is a modification of the standard WKB quantization for ob- 
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taining the quasiclassical eigenvalues of nonrelativistic Hamiltonians. It was demonstrated 
that the leading-order SWKB quantization condition in each and every case reproduces 
the exact energy eigenvalues for a class of solvable potentials. For these models, solutions 
can be written in terms of elementary functions. 

Successes of the SWKB quantization rule have revived interest in the original WKB 
quantization condition. In several common applications the method gives very accurate 
results. Proofs of varying degrees of rigor have been advanced that demonstrate the 
exactness of the standard WKB quantization condition The existing proofs 



of exactness of the WKB approximation are not entirely rigorous since the correction 
terms are only asymptotically valid, i.e., as h —>■ flO| . Furthermore, in the cases when a 
modified WKB integral gives the exact eigenvalues, it is not even clear which "correction" 
must be shown to be zero. 

The standard lowest-order WKB prescription reproduces the exact energy levels for 
the harmonic oscillator in the Cartesian coordinates x, y, and z. But just this problem is 
correctly formulated in the framework of the quasiclassical approach; in the Cartesian co- 
ordinates, the Schrodinger's equation has the required canonical form and the generalized 
moments for each degree of freedom coincide with the corresponding classic moments. As 
for other coordinate systems, for example spherical, the WKB method does not reproduce 
the exact energy levels unless one supplements it with Langer-like correction terms. 

For the central potential V(r), the Schrodinger's equation can be written in the spher- 
ical coordinates as 
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r 2 sin 2 6 dip 2 

2m[E - V(r)]if)(f). 



The standard solution of this equation is the following. If one substitutes ip(r) = [U(r)/r] 
[©(0)/Vsin 9]$(ip) into Eq. ([[]), one obtains (after separation) the following three reduced 
one-dimensional equations 



tf-f- + 2m(E-V)- , 

ar z ^ 



L 2 
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ti 



+ L 2 + 
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U(r) = 0, 
0(9) = 0, 



h 2 £- 2 + Llj *{<p) = 0. 



(2) 

(3) 
(4) 



Exact solution of Eq. (||) gives, for the squared angular momentum L 2 , L 2 = 1(1 + l)h 2 . 
Application of the leading-order WKB quantization condition |] , 



/ Jp 2 (x, E)dx = nh (n + - 
J xi V 2 



71=0,1,2,..., 



(5) 
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to the radial Eq. (Q) does not reproduce the exact energy spectrum {here in ([|) X\ and 
%2 are the classic turning points and p 2 (x, E) = 2m[E — V(x)]}. The problem comes from 
the form of the centrifugal term, 1(1 + l)h 2 /r 2 . 

To overcome this problem in particular case of the Coulomb potential, a special tech- 
niques has been developed. In order for the first-order WKB approximation to give the 
exact eigenvalues, the quantity 1(1 + 1) in Eq. @ must be replaced by (/ + |) 2 fill]. The 



reason for this modification (for the special case of the Coulomb potential) was pointed 
out by Langer (1937) |]TIJ from the Langer transformation r = e x , U(r) = e x ^ 2 X(x). 
However, for other spherically symmetric potentials, in order to obtain the appropriate 
Langer-like correction terms, another special transformation of the wave function (w.f.) 
and its arguments is required. 

There are several other related problems in the semiclassical consideration of the radial 
Schrddinger equation (|2J). (i) The WKB solution of the radial equation is irregular at 

r — ► 0, i.e. R WKB (r) oc r x /y/r, X = J 1(1 + 1), whereas the exact solution in this limit 



is R(r) oc r l . (ii) Equation (|2]) has no the centrifugal term when I = 0, i.e. the radial 
problem has only one turning point and one can not use the WKB quantization condition 
(13) derived for two-turning-point problems. However, solving the equation for / > by 
known exact methods one obtains energy eigenvalues for all I. (iii) The WKB solution of 
equation (|3[) has analogous to the radial one, incorrect behavior at 9 — > 0: Q WKB (9) oc 9 11 , 
a 2 = m 2 — h 2 /A, while the exact regular solution in this limit is 0["(#) oc #' m L Angular 
eigenfunction Yqq(9, (p) = const, i.e. no nontrivial solution exists. 

As practical use shows the standard leading-order WKB approximation always repro- 
duces the exact spectrum for the solvable spherically symmetric potentials V(r) if the 
centrifugal term in the radial Schrodinger's equation has the form (I + ^) 2 h 2 /r 2 . As will 
be shown below the centrifugal term of such a form can be obtained from the WKB solu- 
tion of equation (|l|) if separation of this three-dimensional equation has performed with 
the help of the correspondence principle. 



3. Separation of the Schrodinger's equation 



There are two essential features of the WKB method. First, the method was developed 
to solve the Schrodinger's equation in canonical form (without first derivatives). Second, 
in the quasiclassical method, the classic quantities such as classic momentum, classic 
action, phase, etc., are used. (For example, in the WKB quantization condition, the 
classic generalized momentum in the phase-space integral is used). For the harmonic 
oscillator in the Cartesian coordinates, the generalized moments in the original equation 
and moments obtained after separation coincide with the corresponding classic moments. 
Just for this problem, the standard WKB method (in one and multi-dimensional cases) 
reproduces the exact energy levels without any additional correction terms. 

The generalized moments in Eqs. (@)-@ obtained from separation of Eq. ([!]) are 
different from the corresponding classic moments. As a result, the WKB method does 
not reproduce the exact energy levels for the spherically symmetric potentials (unless one 
supplements it with Langer-like correction terms). The reason is the form of the squared 
angular momentum, L 2 = 1(1 + l)h 2 , which is obtained from solution of the equation 
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(3). In the WKB method, in order to reproduce the exact energy spectrum, the term 
(I + \) 2 Ti 2 should be used in the centrifugal term. This term, M 2 = L 2 + h/A, is in the 
angular equation (|3|), but is not in the radial equation @. 

Let us show that the term (I + \) 2 Ti 2 can be obtained from the quasiclassical solution 
of the reduced Schrodinger's equation. For this, exclude in Eq. (|l|) the first derivatives 
that can be easily done with the help of the following operator identity: 
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Then, after dividing by ip(r) = i?(r) 0(0) $(<£>), where R(r) = rR(r), 0(0) 
&((p) = $ (</>), we obtain the equation 
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we can write Eq. (ffi) in the form of the classic Hamilton- Jacobi equation, 
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where So — S (r , E) is the classic action of the system. 

Using the correspondence principle, we see, from Eq. (|TT|), that Eqs. (|8|)-(|10|) are the 
squared generalized moments expressed via the quantum-mechanical quantities. Now, 
let us separate equation (|TTD. Then, taking into account Eqs. (jSP-fllPp, we obtain the 
following system of the second-order differential equations in canonical form 
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where M 2 and M 2 are the constants of separation and, at the same time, integrals of 
motion. 

Equations (@)-(0) have the quantum- mechanical form ftp = fip, where / is the phys- 
ical quantity (the squared generalized momentum) and / is the corresponding operator. 
We see that the term h 2 /4 in the left-hand side of the equations is related to the squared 
angular momentum operator. This term disappears in the leading h approximation H 



and the equations (12|)-(14) can be written in the general form as 




= {\ 2 -U{qmq). (15) 

The squared generalized moments in the right-hand sides are the same as the classic ones. 
The use of these moments in the WKB quantization condition and WKB solution yields 
the exact energy spectra for the central-field potentials and does not result in the diffi- 
culties of the WKB method mentioned above. 

4. Solution of the Schrodinger's equation 

In this section we consider quasiclassical solution of the Schrodinger's equation for the 
spherically symmetric potentials. In order for the approach we consider here to be self- 
consistent we have to solve each equation obtained after separation by the same method, 
i.e. the WKB method. The quasiclassical method is general enough and the WKB for- 
mulas can be written differently, i.e. on the real axis [|IJ and in the complex plane 0. 
Most general form of the WKB solution and quantization condition can be written in the 
complex plane. 

The WKB quantization in the complex plane. The WKB method is usually used 
to solve one- dimensional two turning point problems. Within the framework of the WKB 
method the solvable potentials mean those potentials for which the eigenvalue problem 
has two turning points. However the WKB method can be used to solve problems with 
more then two turning points. In this case formulation in the complex plane is the most 
appropriate. 

Consider Eqs. (P^ )-(|14T) in the framework of the quasiclassical method. Solution of 
each of these equations [in the general form Eq. (0)] we search in the form 



ij)x(z) = A exp 



k s{zA) 



(16) 



where A is the arbitrary constant. The function S(z, A) is written as the expansion in pow- 
ers of h, S(z, A) = Sq(z, A) + hSi(z, A) + \h 2 S2{z 1 A) + . . .. In the leading h approximation 
the WKB solution of Eqs. flT2|)-(|l4]) can be written in the form 



j, WKB {z) = A exp 



(17) 



In quantum mechanics, quantum numbers are determined as number of zeroes of the 
w.f. in the physical region. In the complex plane, the number of zeroes N of a function 
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y(z) inside the contour C is defined by the argument's principle |T5|, [16). For the w.f. 
ipx(z), according to this principle we have 



c i/;\(z) 



dz = 2niN, 



(18) 



where ip'\(z) is the derivative of the function ip\(z) over the variable z [see Ref. [|T7|] for 



more information about the condition (Tg)]. Contour C is chosen such that it includes cuts 
(therefore, zeroes of the w.f.) between the turning points where p 2 (z, A) = A 2 — U(z) > 0. 
Substitution of Eq. 



into (|I~8"D results in the quantization condition 



2TihN. 



(19) 



In the case, when p(z, A) is a smooth function of the spatial variable and the equation 
A 2 — U(z) = has two roots (turning points), the quantization condition ( |T9"D takes the 
form 



j ^p 2 (z,X)dz = 2nh (n+^ 



In particular, for p 2 (z, A) = A 2 , the quantization condition is 



^p 2 (z,X)dz 



2<k%N. 



(20) 



(21) 



In the next section we solve the three-dimensional Schrodinger equation for several spher- 
ically symmetric potentials by the method under consideration. 



A. The angular momentum eigenvalues. Equations ( |i~3l) and (|i"4l) determine the 
squared angular momentum eigenvalues, M 2 , and its projection, M z , respectively. The 



quantization condition (|2TD appropriate to the angular equation (JlJ) 



<j> M z d(f = 2-Khm, 

^ives M z = km, m = 0, 1, 2, ... The corresponding quasiclassical solution is 



$ m (^) = C , 1 e im ^ + C 2 e- 



irrup 



[22) 



(23) 



where C\ and C 2 are the arbitrary constants. 

The quantization condition (|^) appropriate to Eq. flilf ) is 




Ml 



dO = 2nh [n e + 



n 6 = 0,1,2, 



(24) 



sin 2 fl"~ "" V" ' 2, 

To calculate the integral (2^) (as other hereafter) we use the method of stereographic 
projection. This means that, instead of integration about a contour C enclosing the 
classical turning points, we exclude the singularities outside the contour C, i.e., at 9 = 
and oo in this particular case. Excluding these infinities we have, for the integral (|24|), 
I — Iq + Too. Integral Iq = —27rM z , and 1^ is calculated with the help of the replacement 
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e %e that gives I 00 = 2tx^ 



2nM. Therefore, I = 2ir(M — M z ) and we obtain, for 



the squared angular momentum eigenvalues, 



M 



1 \ 2 



(25) 



where I = n e + m. Thus the quasiclassical solution of the Schrodinger's equation results 
in the squared angular momentum eigenvalues (|23|). This means the centrifugal term in 
the radial Eq. ( |12"D has the form (/ + \) 2 fi 2 jr 2 for all spherically symmetric potentials. 

As known the WKB solution Q WKB (9) of the equation (Q) has incorrect asymptotes 
at 9 — > and n. At the same time, the WKB solution of Eq. (13), which corresponds 
to the eigenvalues (pS|), has the correct asymptotic behavior at these points for all I. So 
far, as the generalized momentum p(9) ~ at 9 — » 0, this gives, for the WKB solution 
in the representation of the wave function tp(r), ©["($) = Q WKB (9)/y/sin 9 oc # H which 
corresponds to the behavior of the known exact solution Yi m (9, (p) at 9 — > 0. 

In the classically allowed region, where p 2 (9, M) = M 2 — M 2 / sin 2 9 > 0, the leading- 
order WKB solution of Eq. (|I3]) is 



e 



WKB 



(0) 



B 



p(9,M) 



: COS 



p(9,M)d9- 



7T 



(26) 



The normalized quasiclassical solution far from the turning points, where p(9, M) 
(I + |)?t, can be written in elementary functions as 



(0) 



21 + 1 



\ 7r(Z - m + \ 



■ cos 



'4 



-(/-m) 



(27) 



where we have took into account that the phase-space integral at the classic turning 
point 9i is x(#i) = — f (ne + \) and x(^) = \{n e + |) at 6 1 = We see that the 
eigenfunctions (|27D are either symmetric or antisymmetric. The corresponding WKB 
solution, Y^ KB (9, <p) = 6J n (^)$ m (v9), where the normalized eigenfunction $ m (v?) ~ 1 



■2jt 



e ±imip , in the representation of the w.f. ip(r) is 



Y, 



WKB 
Im 



1 



m + 2 



COS 



7T 



(1 



m 



(28) 



Remind some results concerning the semiclassical approach in quantum mechanics. 
The general form of the semiclassical description of quantum-mechanical systems have 
been considered in Ref. | 18| . It was shown that the semiclassical description resulting 
from Weyl's association of operators to functions is identical with the quantum description 
and no information need to be lost in going from one to the another. What is more "the 



IS 



semiclassical description is more general than quantum mechanical description..." 
The semiclassical approach merely becomes a different representation of the same algebra 
as that of the quantum mechanical system, and then the expectation values, dispersions, 
and dynamics of both become identical. 
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One of the fundamental features of quantum mechanical systems is nonzero minimal 
energy which corresponds to quantum oscillations. The corresponding w.f. has no zeroes 
in the physical region. Typical example is the harmonic oscillator. 

Eigenvalues of the one-dimensional harmonic oscillator are E n = fuo{n + |), i.e. the 
energy of zeroth oscillations E = ^hu). In three-dimensional case, in the Cartesian 
coordinates, the eigenvalues of the oscillator are E n = Tiuj{n x + n y + n z + |) [[HJ, i.e. each 
degree of freedom contributes to the energy of the ground state, E = E 0tX + Eq )V + E 0tZ = 
^hu. Energy of the ground state should not depend on coordinate system. This means 
that, in the spherical coordinates, each degree of freedom (radial and angular) should 
contribute to the energy of zeroth oscillations. In many applications and physical models 
a nonzero minimal angular momentum M is introduced (phenomenologically) in order to 
obtain physically meaningful result (see, for instance, Ref. |2(|). However, the existence 
of M follows from the quasiclassical solution of Eq. (|T3"D 0, jZTJ . 

Consider the WKB eigenfunction (|28|) for the ground state. Setting in (|28|) m = and 
/ = 0, we obtain the nontrivial solution in the form of a standing half- wave [remind that 
the spherical function Y 0Q (9, <p) = const], 

Y J KB (9,<f) = -cos 9 -. (29) 

7T Z 

The corresponding eigenvalue is 

M = |. (30) 

The eigenvalue (p0|) contributes to the energy of zeroth oscillations. This means that (p9|) 
can be considered as solution, which describes the quantum fluctuations of the angular 
momentum. Note, that the eigenfunction of the ground state, Y^ KB (9, ip), is symmetric. 



Below, we solve the radial equation (|12[ ) for some spherically symmetric potentials and 



show the contribution of the eigenvalue M to the energy of the ground state. 

B. The Coulomb problem V(r) = —-. The WKB quantization condition (19) 
appropriate to the radial equation (0) with the Coulomb potential does not reproduce the 
exact energy levels unless one supplements it with Langer-like correction terms. Another 
problem is that the radial Schrodinger's equation (Q) has no the centrifugal term at / = 
and one can not use the WKB quantization condition (derived for two-turning-point 
problems) to calculate the energy of the ground state directly from this equation. We 
do not run into such a problem in case of Eqs. ([12])- (|T^). As follows from the above 
consideration, the centrifugal term (/ + ^) 2 h 2 /r 2 is the same for all spherically symmetric 
potentials and the WKB method reproduces the exact energy spectrum for all / and n r . 

For the radial equation (|I^) with the Coulomb potential, the WKB quantization con- 
dition (EOl) is 



c 



2ma M 2 , / 1\ 

2mE + -dr = 2-kU (n r + - , (31) 



\ Lmni -\ —ar = znn \ n r -\ — 

\ r r 2 \ 2, 

where the integral is taken about a contour C inclosing the turning points r± and r 2 . Using 
the method of stereographic projection, we should exclude the singularities outside the 
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contour C, i.e. at r = and oo. Excluding these infinities we have, for the integral (|3l|). 



I = Iq + loo, where Iq = 2ni\J — M 2 = —2nM and 1^ = 2niam / y2mE . The sequential 
simple calculations result in the exact energy spectrum 



a 2 m 



En ~ 2[{n r + \)h + M} 2 ' (32) 

For the energy of zeroth oscillations we have, from Eq. (|32|), E = — ^a 2 m(^ + M ) -2 , 
that apparently shows the contribution of the quantum fluctuations of the angular momen- 



tum (see Eq. fl3"0|)) into the energy of the ground state E pifl . The radial quasiclassical 
eigenfunctions, R^ KB (r), inside the classical region [ri,r2] far from the turning points r\ 
and T2 are written in elementary functions in the form of a standing wave H, 

R n (r) = Acos (\p n r + ^ri r ) , (33) 



X 2 

where we have took into account that the phase-space integral ([$!]) at the classic turning 
point r\ is x( r i) = —j( n r + \)- Here A is the normalization constant and p n is the eigen- 



momentum expressed via the energy eigenvalue E n , p n = J2m\E n \. The eigenfunctions 



fl33|) are either symmetric or antisymmetric. 



C. The three-dimensional harmonic oscillator V(r) = \mjj 1 r 1 . The three- 
dimensional harmonic oscillator is another classic example of the exactly solvable problems 
in quantum mechanics. The problem has 4 turning points, r±, r2, r%, and r^, but only two 
of them, r 3 and r 4 , lie in the physical region r > 0. 

The problem is usually solved with the help of the replacement x = r 2 which reduces 
the problem to the 2-turning-point (2TP) one. But this problem can be solved as the 
4TP problem in the complex plane. Because of importance of the oscillator potential in 
many applications and with the purpose of further development of the WKB method, we 
shall solve the problem by two methods, on the real axis as 2TP problem and then in the 
complex plane as 4TP problem. 

Consider first the physical region r > 0, where the problem has two turning points. 
The leading-order WKB quantization condition (p0|) then is 



/• r 4 M 2 ( 1\ 

J \ 2mE — {mur) 2 — dr = nh yn r + -J , (34) 



where n r is the number of zeroes of the w.f. between the classic turning points r 3 and 
r 4 . Integral (|34"D is reduced to the above case of the Coulomb potential with the help of 
the replacement z = r 2 . Integration result is I = tt(E/uj — M)/2 and we obtain, for the 
energy eigenvalues, 



E n = UJ 



27) [n r + -)+ M 



(35) 



So far, as M = (I + ^)h, we obtain the exact energy spectrum for the isotropic oscillator. 
Energy of the ground state is E = uj(h + M ), where M is the contribution of quantum 
fluctuations of the angular momentum. 
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Emphasize the following in this solution. The 4TP problem has been solved as the 2TP 
problem; we have applied the 2TP quantization condition (p0| ) to the 4TP problem that is 
not quite correct. We have obtained the correct result because the potential is symmetric 
and the replacement x = r 2 reduces the problem to the 2TP problem, i.e. "reflects" the 
negative region r < (and zeroes of the w.f.) into the positive region. A more correct 
approach to solve the problem would be a 4TP quantization condition. Fortunately, the 
WKB method in the complex plane allows to solve this problem as the 4TP problem. 

In the complex plane, the problem has two cuts, between turning points r 1; r 2 and r 3 , 
r4. To apply residue theory for the phase space integral we need to take into account all 
zeroes of the w.f. in the complex plane, i.e. the contour C has to include both cuts. The 
quantization condition ([19]) in this case takes the form 



c 



p(r, E) + i 



Ci 



.hp'jr, E) 



dr + f 

Jc 2 



.h p'(r,E) 



hp'{r,E) 
2 p(r, E) 

dr = 



dr = 
2-kUN, 



(36) 



where p 2 (r, E) = 2mE — (mur) 2 — M 2 /r 2 , and C\ and C% are the contours about the cuts 
at r < and r > 0, respectively. The number N = n r<0 + n r>0 , where n r<0 and n r>0 
are the numbers of zeroes of the w.f. at r < and r > 0, respectively. For the harmonic 
oscillator, because of symmetricity of the potential we have n r< o = n r> o = n r , i.e. the 
total number of zeroes is N = 2n r . 

Therefore, the quantization condition (Sol) for the 4TP problem takes the form, 



n r + 



(37) 



f p(r, E)dr = f p(r, E)dr + f p(r, E)dr = 2-Khk 

Jc JCi JC-2 

where k = 2 is the number of cuts. We can write the 4TP quantization condition in this 
form because the effective potential is infinite at r = 0. In case if the potential is finite in 
the whole region, the quantization condition will be more complicate . 

The condition ( |3"TD is in agreement with the Maslov's theory. This means that the 
right-hand side of the equation (|37|) can be written in the form 



2nhk I n 



(3? 



where \x = 2k is the Maslov's index, i.e. number of reflections of the w.f. on the walls of 
the potential. 

In the general case of the potential which is infinite between cuts, the 2k turning point 
quantization condition is 



I p(z,E)dz = 2nhY^ 

^ C 8=1 



n 



2nh[N + 



(39) 



where N = krii is the total number of zeroes of the w.f. on the k cuts. On the real axis, 
for the 2k TP problem, the quantization condition (39) has the form 
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E / ^P 2 (z,E)dz = nh(N + n. (40) 

Because the harmonic oscillator potential is symmetric, integrals in Eq. ( |3"TD are iden- 
tical, i.e. the quantization condition (|37l) is equivalent to the 2TP quantization condition 
fl2"0"|). The phase-space integral can be easily calculated in the complex plane. For this 
we have to exclude the singularities at r = and oo outside the contour C . Excluding 
these infinities we have, for the integral (|39|) with the isotropic potential, / = I + 1^, 
where Iq = —2ttM and integral 1^ is calculated with the help of the replacement r = l/z, 
Joe = 2ttE/u, i.e. we again obtain the exact result (^) for E n . 

Consider Eq. fl35|) at n r = and / = 0, i.e. the energy of the ground state. We have 
E = u>(H + M ), where M = Ti/2 is the contribution of the quantum fluctuations of 
the angular momentum into the energy of the ground state E . The radial quasiclassical 
eigenfunctions, R^ KB (r), in the region of the classical motion far from the turning points 
are written analogously to the above case in the form of a standing wave [see Eq. (^)] . 

D. The Hulthen potential V(r) = -V e- r/r °/(l - e~ r/ro ). The Hulthen potential 
is of a special interest in atomic and molecular physics. The potential is known as an 
"insoluble" by the standard WKB method potentials, unless one supplements it with 
Langer-like corrections. The radial problem for this potential is usually considered at 
I = 0. However, in the approach under consideration, the quasiclassical method results in 
the nonzero centrifugal term at / = and allows to obtain the analytic result for all I. 

The leading-order WKB quantization condition ( p0|) for the Hulthen potential is 



/ e -r/r \ M 2 / 1\ 

f m [ E + Vo x-e-r/ro ) - ^ dr = M + 2 ) ■ (41) 

In the region r > 0, this problem has two turning points r\ and r 2 . The phase-space 
integral Q4"l"| ) is calculated analogously to the above case. Introducing the new variable 
p = r/r , we calculate the contour integral in the complex plane, where the contour 
C encloses the classical turning points pi and p2- Using the method of stereographic 
projection, we should exclude the infinities at r = and oo outside the contour C. 
Excluding these infinities we have, for the integral (f4lD , I = I + 1^, where Io = —2tiM 
and Too is calculated with the help of the replacement z = e p — 1 0, 



2,nr? ) [E- \ - — ] - —dp = (42) 



-2-kM + 27rr V z 2m 



-V^E + J-E + V Q 



Substituting the integration result into Eq. (|4l"|), we immediately get the exact energy 
spectrum 

(43) 
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where N = (n r + + M is the principal quantum number. Setting in fl43|) M = 0, we 
arrive at the energy eigenvalues obtained from known exact solution of the Schrodinger's 
equation at / = 0. However, in our case M min = M = h/2 at I — and the principal 
quantum number is N = (n r + \)h + M . As in the previous examples, this apparently 
shows the contribution of the quantum fluctuations of the angular momentum into the 
energy of the ground state, E . 

E. The Morse potential V(r) = V r [e~ 2a(r/ro-1) -2 e - a ^ r °-% The Morse potential 
is usually considered as one- dimensional problem at I = 0. In this case the problem has 
two turning points (note that the left turning point, r\, is negative) and can be solved 
exactly. In the general case, for I > 0, we have an "insoluble" 4TP problem. 

For this potential, let us consider, first, the radial Schrodinger equation (Q), which 
does not contain the centrifugal term at I = 0, 

(~ l7i £^) U ^ = 2m [ E ~ V e- 2a{r - ro)/ro + 2V e- a{r - ro)/ro ] U(r). (44) 
The first-order WKB quantization condition (||) appropriate to this equation is 

\j2m\E — y e- 2a (r-n))/ro + 2V Q e-< r - r ^/ r ^]dr = nh (n r + -) . (45) 

ri V 2/ 

Introducing a variable X = e ~«(»-n))/ro ) we 

reduce the phase-space integral to the well 
known one. Sequential simple calculations result in the exact energy eigenvalues 



E n = -Vq 



1 - 



ah(n r + \] 



(46) 



r y/2mV 

Now, let us deal with Eq. (^) for this potential, which contains the non-vanishing 
centrifugal term, ^. 2 /4r 2 , at I = 0. In this case we have an "insoluble" 4TP problem. In 
the complex plane, the problem has two cuts (k = 2), at r < and r > 0, therefore, we 
apply the 4TP quantization condition (|37|), 



c 



\ 



2m [E - l/ e- 2 °( r - r °)/ r( > + 2Voe- a ( r - r °V r °] 



M 2 



dr = 47rh { n r + 



(47) 



where the contour C encloses the two cuts, but does not enclose the point r = 0. To 
calculate this integral introduce the variable p = r/r . Using the method of stereographic 
projection, we should exclude the singularities outside the contour C, i.e. at r = and 
oo. Excluding these infinities we have, for the integral ((47|) ||, 



-2nM - ( V-2mE 

a 



2mVn 



(48 



and for the energy eigenvalues this gives 



1 — a 



2h{n r + \) + M 
r y/2mV 



(49) 
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Setting in fl49|) I = 0, we obtain, 



a[2h(n 



E n = -V Q !_ "^ + 2 + M, . (5Q) 

Equation ( |50l) for is different from the expression fl46|) obtained from solution of Eq. 
(H) for the Morse potential at I — 0. This difference is caused by the nonzero centrifugal 
term h 2 /Ar 2 in the radial equation (|l|) at I — 0. Thus we obtain two results for the Morse 
potential by the WKB method: the known exact eigenvalues (fllf ) obtained from solution 
of Eq. (H) at / = and another result (|49| ) obtained from solution of Eq. (|T2|) for all /. 



1 2 



F. The potential V(r) = kr + |u; 2 r 2 . This linear plus isotropic potential is one of 
the interest in particle physics. The potential has four turning points and can not be 
reduced to a hypergeometric function by a suitable transformation. This multi-turning- 
point problem is "insoluble", also, by the standard WKB method. However, the WKB 
method in the complex plane allows easily to solve this 4TP problem [see, for instance, 
Refs. H, H]. 

The problem has two cuts (k = 2) in the complex plane between turning points r 1; r 2 
and r 3 , r 4 . Quantization condition (|39|), for this problem, is 



c 



\ 



2mE — 2mkr 



M 2 



\muor} 



dr = Anh { n r + 



(51) 



where the contour C includes both cuts, but includes no the point r = 0. To calculate the 
integral (^T|) we exclude the infinities at r = and oo outside the contour C. Excluding 
these infinities we have, for the integral (pT|), 



2tt 



or, for the energy eigenvalues, this gives 




M 



Anh n r + 



(52) 



E„ 



2h[n r + - ) + M 




(53) 



Energy eigenvalues ( [5*3"D are similar to the harmonic oscillator ones but shifted by the 
constant value —-^{k/u) 2 . Putting in Eq. ( |53"D k = 0, we arrive to the eigenvalues (|35|) 
for the isotropic oscillator. 

Analogously one can obtain energy eigenvalues for other spherically symmetric poten- 
tials. The standard leading-order WKB approximation appropriate to the wave equation 
( |T2"D yields the exact energy eigenvalues for known solvable potentials and "insoluble" 
ones with more than two turning points. This is possible because the centrifugal term 
of the required form, (/ + \) 2 fi /r 2 , has obtained in a natural way from solution of the 
angular equation (|T3| ) with the use of the same WKB method; this term is the same for 
all central-field potentials. 
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5. Conclusion 

Conventional approach to solve the Schrodinger's equation is to reduce the original 
equation to a hypergeometric form or some special function by a suitable transformation. 
In each and every case, one needs to find, first, a special transformation for the wave 
function and its arguments to reduce the original equation to some known equation. 
There is another way to solve the Schrodinger's equation which is simple, general for all 
types of problems in quantum mechanics, and a very efficient to solve not only two- but 
multi-turning point problems. 

Almost together with quantum mechanics an appropriate method to solve the wave 
equation has been developed known mainly as the WKB approximation. This method is 
general for all types of problems in quantum mechanics, simple from the physical point 
of view, and its correct application results in the exact energy eigenvalues for all solvable 
potentials. 

In spite of long history no any strict rules concerning the application of the WKB 
method to multi-dimensional problems in quantum mechanics have been formulated. 
Meanwhile, this topic is closely related to the problem of exactness of the WKB method. 
The exactness of the method has proven in the literature for many potentials with the help 
of specially developed techniques, or improvements, or modifications of the quasiclassical 
method on the real axis and in the complex plane. In this work we have fulfilled the qua- 
siclassical analysis of the three-dimensional Schrodinger's equation. The original equation 
has been reduced to the form of the classic Hamilton- Jacobi equation without first deriva- 
tives. Separation of the equation has been performed with the help of the correspondence 
principle between classic and quantum-mechanical quantities. As a result of the separa- 
tion, we have obtained the system of reduced second-order differential equations. Each 
of these equations has the correct quantum-mechanical form, p 2 ip(q) = p 2 {q)ip{q), and 
solved by the WKB method. We have stressed that the squared generalized moments, 
p 2 (q), obtained after separation should coincide with the corresponding classic moments. 
This means that the problem under consideration should correspond to a concrete classic 
problem. 

We have shown that the Langer replacement 1(1 + 1) — > (I + \) 2 needed to repro- 
duce the exact energy spectrum for the spherically symmetric potentials by the WKB 
method requires the modification of the squared angular momentum in the quasiclassical 
region. The squared angular momentum eigenvalues, M 2 = (I + \) 2 h 2 \ have obtained in 
our approach from solution of the angular wave equation in the framework of the same 
quasiclassical method. As a result, the centrifugal term has the form (1 + \) 2 fi 2 /r 2 for any 
spherically symmetric potential V(r). 

The quasiclassical solution contains a more detail information in comparison with 
known exact solution. One of consequences of the WKB solution of the Schrodinger's 
equation is the existence of a nontrivial angular eigenfunction of the type of a standing 
half- wave for the angular quantum number / = 0. This solution has treated as one 
which describes the quantum fluctuations of the angular momentum with the eigenvalue 
M = h/2. We have shown that the quantum fluctuations of the angular momentum 
contribute to the energy of the ground state, E . 

To demonstrate efficiency of the quasiclassical method, we have solved the three- 
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dimensional Schrodinger's equation for some central-field potentials. The quasiclassical 
method successfully reproduces the exact energy spectrum not only for solvable spheri- 
cally symmetric potentials but, also, for "insoluble" potentials with more than two turning 
points. The quasiclassical eigenfunctions for the discrete spectrum have written in ele- 
mentary functions in the form of a standing wave. 

The remarkable features of the quasiclassical method incline us to treat the leading- 
order WKB approximation as a special (asymptotic) exact method to solve the Schrodinger 
equation. In the quasiclassical approach we use the same technique for all types of prob- 
lems. The same simple rules formulated for two-turning problems work for many turning 
point problems, as well. In this sense, the quasiclassical method is a more general in 
comparison with traditional one with the use of techniques of the special functions. 
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